Chapter 1 



PHYSICS OF NEUTRON STAR KICKS 

Dong Lai 

Center for Radiophysics and Space Research, Cornell University, Ithaca, NY 14853, USA 
dong@spacenet.tn.cornell.edu 

Abstract It is no longer necessary to "sell" the idea of pulsar kicks, the notion that neutron 
stars receive a large velocity (a few hundred to a thousand km s _1 ) at birth. 
However, the origin of the kicks remains mysterious. We review the physics 
of different kick mechanisms, including hydrodynamically driven, neutrino and 
magnetically driven kicks. 

1. INTRODUCTION 

It has long been recognized that neutron stars have space velocities much 
greater than their progenitors'. Recent studies of pulsar proper motion give 
200 — 500 km s _1 as the mean 3D velocity of NSs at birth (e.g., Lyne and 
Lorimer 1994; Lorimer et al. 1997; Hansen & Phinney 1997; Cordes & Chernoff 
1998), with possibly a significant population having velocities greater than 
700 km s _1 . Direct evidence for pulsar velocities >1000 km s" 1 comes 
from observations of the bow shock produced by the Guitar Nebula pulsar 
(B2224+65) in the interstellar medium (Cordes et al. 1993) and studies of 
pulsar-supernova remnant associations (e.g., Frail et al. 1994; Kaspi 1999). A 
natural explanation for such high velocities is that supernova explosions are 
asymmetric, and provide kicks to nascent neutron stars. Support for supernova 
kicks has come from the detection of geodetic precession in binary pulsar PSR 
1913+16 (Cordes et al. 1990; Kramer 1998; Wex et al. 1999) and the orbital 
plane precession in PSR J0045-7319/B-star binary (Kaspi et al. 1996) and 
its fast orbital decay (which indicates retrograde rotation of the B star with 
respect to the orbit; Lai 1996); These results demonstrate that binary break- 
up (as originally suggested by Gott et al. 1970; see Iben & Tutukov 1996) 
can not be solely responsible for the observed pulsar velocities, and that natal 
kicks are required. Evolutionary studies of neutron star binary population also 
imply the existence of pulsar kicks (e.g., Fryer et al. 1998). Finally, there are 
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many direct observations of nearby supernovae (e.g., Wang et al. 1999) and 
supernova remnants which support the notion that supernova explosions are 
not spherically symmetric. 

Despite decades of theoretical investigations, our understanding of the phys- 
ical mechanisms of core-collapse supernovae remains significantly incomplete. 
While there is a consensus that neutrino heating of the stalled shock (0.1 — 1 s 
after bounce) plays an important role in driving the explosion, it is unclear 
whether the heating is sufficient to produce the observed supernova energetics; 
It is also unclear whether any convective motion or hydrodynamical instability 
is central to the explosion mechanism (e.g., Herant et al. 1994; Burrows et 
al. 1995; Janka & Miiller 1996; Mezzacappa et al. 1998). The prevalence of 
neutron star kicks poses a significant mystery, and indicates that large-scale, 
global deviation from spherical symmetry is am important ingredient in any 
successful theory of core-collapse supernovae. 

In this paper, we concentrate on different classes of physical mechanisms for 
generating neutron star kicks (§§2-4), and then briefly discuss the astrophysi- 
cal/observational implications (§5). 

2. HYDRODYNAMICALLY DRIVEN KICKS 

The collapsed stellar core and its surrounding mantle are susceptible to a 
variety of hydrodynamical (convective) instabilities (e.g., Herant et al. 1994; 
Burrows et al. 1995; Janka & Miiller 1996; Keil et al. 1996). It is natural to ex- 
pect that the asymmetries in the density, temperature and velocity distributions 
associated with the instabilities can lead to asymmetric matter ejection and/or 
asymmetric neutrino emission. Numerical simulations, however, indicate that 
the local, post-collapse instabilities are not adequate to account for kick veloci- 
ties >100 km s" 1 (Janka & Miiller 1994; Burrows & Hayes 1996; Janka 1998) 
— These simulations were done in 2D, and it is expected that the flow will be 
smoother on large scale in 3D simulations, and the resulting kick velocity will 
be even smaller. 

Global asymmetric perturbations in presupernova cores are required to pro- 
duce the observed kicks hydrodynamically (Goldreich et al. 1996). Numerical 
simulations by Burrows & Hayes (1996) demonstrate that if the precollapse 
core is mildly asymmetric, the newly formed neutron star can receive a kick 
velocity comparable to the observed values. (In one simulation, the density of 
the collapsing core exterior to O.9M and within 20° of the pole is artificially 
reduced by 20%, and the resulting kick is about 500 km s _1 .) Asymmetric 
motion of the exploding material (since the shock tends to propagate more 
"easily" through the low-density region) dominates the kick, although there is 
also contribution (about 10 — 20%) from asymmetric neutrino emission. The 
magnitude of kick velocity is proportional to the degree of initial asymmetry in 
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Figure 1.1 Propagation diagram computed for a 15M Q presupernova model of Weaver and 
Woosley (1993). The solid curve shows iV 2 , where N is the Brunt-Vaisala frequency; the 
dashed curves show L 2 , where Li is the acoustic cutoff frequency, with 2 = 1,2, 3. The spikes 
in N 2 result from discontinuities in entropy and composition. The iron core boundary is located 
at 1.3Mq, the mass-cut at 1.42M . Convective regions correspond to N = 0. Gravity modes 
(with mode frequency ui) propagate in regions where lo < N and w < Li, while pressure modes 
propagate in regions where u> > N and w > L t . Note that a g-mode trapped in the core can 
lose energy by penetrating the evanescent zones and turning into an outgoing acoustic wave (see 
the horizontal line). Also note that g-modes with higher n (the radial order) and I (the angular 
degree) are better trapped in the core than those with lower n and I. 

the imploding core. Thus the important question is: What is the origin of the 
initial asymmetry? 

(i) Presupernova Perturbations 

Goldreich et al. (1996) suggested that overs table g-mode oscillations in the 
presupernova core may provide a natural seed for the initial asymmetry. These 
overstable g-modes arise as follows. A few hours prior to core collapse, a 
massive star (M > 8M ) has gone through a successive stages of nuclear 
burning, and attained a configuration with a degenerate iron core overlaid by 
an "onion skin" mantle of lighter elements. The rapidly growing iron core 
is encased in and fed by shells of burning silicon and oxygen, and the entire 
assemblage is surrounded by a thick convection zone. The nearly isothermal 
core is stably stratified and supports internal gravity waves. These waves cannot 
propagate in the unstably stratified convection zone, hence they are trapped and 
give rise to core g-modes in which the core oscillates with respect to the outer 
parts of the star. The overstability of the g-mode is due to the "e-mechanism" 
with driving provided by temperature sensitive nuclear burning in Si and O 
shells surrounding the core before it implodes. It is simplest to see this by 
considering a I = 1 mode: If we perturb the core to the right, the right-hand- 
side of the shell will be compressed, resulting in an increase in temperature; 
Since the shell nuclear burning rate depends sensitively on temperature (power- 
law index ~ 47 for Si burning and ~ 33 for O burning), the nuclear burning is 
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greatly enhanced; This generates a large local pressure, pushing the core back 
to the left. The result is an oscillating g-mode with increasing amplitude. 

The main damping mechanism comes from the leakage of mode energy. The 
local (WKB) dispersion relation for nonradial waves is 

kl = {^cl)-\^-L^-N% (1.1) 

where k r is the radial wavenumber, L\ = TjUj~+T)c s /r (c s is the sound 
speed) and N are the acoustic cut-off (Lamb) frequency and the Brunt- Vaisala 
frequency, respectively. Since acoustic waves whose frequencies lie above 
the acoustic cutoff can propagate through convective regions, each core g- 
mode will couple to an outgoing acoustic wave, which drains energy from the 
core g-modes (see Fig. 1.1). This leakage of mode energy can be handled 
with an outgoing propagation boundary condition in the mode calculation. 
Also, neutrino cooling tends to damp the mode. Since the nuclear energy 
generation rate depends more sensitively on temperature than pair neutrino 
emission (power law index ~ 9), cooling is never comparable to nuclear heating 
locally. Instead, thermal balance is mediated by the convective transport of 
energy from the shells, where the rate of nuclear energy generation exceeds 
that of neutrino energy emission, to the cooler surroundings where the bulk 
of the neutrino emission takes place. Calculations (based on the 15M Q and 
25M Q presupernova models of Weaver & Woosley 1993) indicate that a large 
number of g-modes are overstable, although for low-order modes (small I and 
n) the results depend sensitively on the detailed structure and burning rates of 
the presupernova models (Lai & Goldreich 2000b, in preparation). 

Our tentative conclusion is that overstable g-modes can potentially grow to 
large amplitudes prior to core implosion, although a complete understanding of 
the global pre-collapse asymmetries is probably out of reach at present, given 
the various uncertainties in the presupernova models (see Bazan & Arnett 1998 
for complications due to convective shell burning in presupernova stars). 

(ii) Amplification of Perturbation During Core Collapse 

Core collapse proceeds in a self-similar fashion, with the inner core shrinking 
subsonically and the outer core falling supersonically at about half free-fall 
speed (Goldreich & Weber 1980; Yahil 1983). The inner core is stable to 
non-radial perturbations because of the significant role played by pressure in 
its subsonic collapse. Pressure is less important in the outer region, making 
it more susceptible to large scale instability. A recent stability analysis of 
Yahil's self-similar collapse solution (which is based on Newtonian theory 
and a poly tropic equation of state P oc p r , with T ~ 1.3) does not reveal any 
unstable global mode before the proto-neutron star forms (Lai 2000). However, 
during the subsequent accretion of the outer core (involving 15% of the core 
mass) and envelope onto the proto-neutron star, nonspherical perturbations 
can grow according to Sp/p oc r -1 / 2 or even 5p/p oc r _1 (Lai & Goldreich 
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2000). The asymmetric density perturbations seeded in the presupernova star, 
especially those in the outer region of the iron core, are therefore amplified (by a 
factor of 5-10) during collapse. The enhanced asymmetric density perturbation 
may lead to asymmetric shock propagation and breakout, which then give rise 
to asymmetry in the explosion and a kick velocity to the neutron star (Goldreich 
et al. 1996; Burrows & Hayes 1996). 

3. NEUTRINO DRIVEN KICKS 

The second class of kick mechanisms rely on asymmetric neutrino emission 
induced by strong magnetic fields. The fractional asymmetry a in the radiated 
neutrino energy required to generate a kick velocity is a = Mv^^c/ 'E tot 
(= 0.028 for t; kick = 1000 km s~\ neutron star mass M = 1.4 Mq and total 
neutrino energy radiated E'tot = 3 x 10 53 erg). 

(i) Effect of Parity Violation 

Because weak interaction is parity violating, the neutrino opacities and 
emissivities in a magnetized nuclear medium depend asymmetrically on the 
directions of neutrino momenta with respect to the magnetic field, and this 
can give rise to asymmetric neutrino emission from the proto-neutron star. 
Chugai (1984) (who gave an incorrect expression for the electron polarization 
in the relativistic, degenerate regime) and Vilenkin (1995) considered neutrino- 
electron scattering, but this is less important than neutrino-nucleon scattering 
in determining neutrino transport in proto-neutron stars. Dorofeev et al. (1985) 
considered neutrino emission by Urea processes, but failed to recognize that in 
the bulk interior of the star the asymmetry in neutrino emission is cancelled by 
that associated with neutrino absorption (Lai & Qian 1998a). 

Horowitz & Li (1998) suggested that large asymmetries in the neutrino flux 
could result from the cumulative effect of multiple scatterings of neutrinos by 
slightly polarized nucleons (see also Janka 1998; Lai & Qian 1998a). However, 
it can be shown that, although the scattering cross-section is asymmetric with 
respect to the magnetic field for individual neutrinos, detailed balance requires 
that there be no cumulative effect associated with multiple scatterings in the 
bulk interior of the star where thermal equilibrium is maintained to a good 
approximation (Arras & Lai 1999a; see also Kusenko et al. 1998). For a given 
neutrino species, there is a drift flux of neutrinos along the magnetic field in 
addition to the usual diffusive flux. This drift flux depends on the deviation of 
the neutrino distribution function from thermal equilibrium. Thus asymmet- 
ric neutrino flux can be generated in the outer region of the proto-neutron star 
(i.e., above the neutrino-matter decoupling layer, but below the neutrinosphere) 
where the neutrino distribution deviates significantly from thermal equilibrium. 
While the drift flux associated with f M 's and zVs is exactly canceled by that 
associated with P„'s and zVs, there is a net drift flux due to z^ e 's and i> e 's. 
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Arras & Lai (1999b) found that the asymmetry parameter for the v e -v e flux is 
dominated for low energy neutrinos (<15 MeV) by the effect of ground (Lan- 
dau) state electrons in the absorption opacity, e a b s — 0.6^15(^/1 MeV) -2 , 
where B\^ = S/(10 15 G), and for high energy neutrinos by nucleon polar- 
ization (~ n m B/T). Averaging over all neutrino species, the total asymme- 
try in neutrino flux is of order a ~ 0.2e a b s , and the resulting kick velocity 
^kick ~ 50 -B15 km s _1 . There is probably a factor of 3 uncertainty in this esti- 
mate. To firm up this estimate requires solving the neutrino transport equations 
in the presence of parity violation for realistic proto-neutron stars. 

(ii) Effect of Asymmetric Field Topology 

A different kick mechanism relies on the asymmetric magnetic field distri- 
bution in proto-neutron stars (see Bisnovatyi-Kogan 1993; However, he con- 
sidered neutron decay, which is not directly relevant for neutrino emission from 
proto-neutron stars). Since the cross section for u e (u e ) absorption on neutrons 
(protons) depends on the local magnetic field strength due to the quantization 
of energy levels for the e~ (e + ) produced in the final state, the local neutrino 
fluxes emerged from different regions of the stellar surface are different. Cal- 
culations indicate that to generate a kick velocity of ~ 300 km s _1 using this 
mechanism alone would require that the difference in the field strengths at the 
two opposite poles of the star be at least 10 16 G (Lai & Qian 1998b). Note that 
unlike the kick due to parity violation [see (i)], this mechanism does not require 
the magnetic field to be ordered, i.e., only the magnitude of the field matters. 

(iii) Exotic Neutrino Physics 

There have also been several interesting ideas on pulsar kicks which rely on 
nonstandard neutrino physics. It was suggested (Kusenko & Segre 1996) that 
asymmetric v T emission could result from the Mikheyev-Smirnov-Wolfenstein 
flavor transformation between v T and v e inside a magnetized proto-neutron 
star because a magnetic field changes the resonance condition for the flavor 
transformation. This mechanism requires neutrino mass of order 100 eV. An- 
other similar idea (Akhmedov et al. 1997) relies on both the neutrino mass and 
the neutrino magnetic moment to facilitate the flavor transformation. More 
detailed analysis of neutrino transport (Janka & Raffelt 1998), however, indi- 
cates that even with favorable neutrino parameters (such as mass and magnetic 
moment) for neutrino oscillation, the induced pulsar kick is much smaller than 
previously estimated (i.e., B » 10 15 G is required to obtain 100 km s _1 kick). 

4. EM RADIATION DRIVEN KICKS 

For completeness, we mention the post-explosion "rocket" effect due to 
electromagnetic (EM) radiation from off-centered magnetic dipole in the pul- 
sar (Harrison & Tademaru 1975). In this mechanism, the neutron star velocity 
comes at the expense of its spin kinetic energy, which is radiated away asym- 
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metrically via EM braking; The kick is attained on the timescale of the initial 
spindown time of the pulsar (i.e., this is not a "natal" kick). The neutron star 
velocity changes according to Mv = eL/c, where L is the EM braking power, 
and e is the asymmetry parameter. Typically, e ~ 0.1(f2s/c)(//<£//z z ), where Q 
is the spin and s is the "off-center" displacement of the dipole; In fact, there is 
theoretical maximum, e = 0.16, achieved for /j,r = 0, /^/Vz = 0. 63 (Qs/c) 
(where /Ur, /x^, fi z are the three cylindrical components of the dipole). The kick 
velocity is along the spin axis, and t^ick ~ 600 (e/0.1)(^o/l kHz) 2 km s _1 
(where uq is the initial spin). Clearly, Even if the neutron star were born with 
maximum rotation rate and e were maintained at near the maximum value, the 
kick velocity would still be at most a few hundred km s _1 . Given that most 
pulsars were born rotating slowly (see Spruit & Phinney 1998), we conclude 
that "EM rocket" cannot be the main mechanism for pulsar kicks. 

5. DISCUSSION 

Statistical studies of pulsar population have revealed no correlation between 
Vkick and magnetic field strength, or correlation between the kick direction and 
the spin axis (e.g., Lorimer et al. 1995; Cordes & Chernoff 1998; Deshpande 
et al. 1999). Given the large systematic uncertainties, the statistical results, 
by themselves, cannot reliably constrain any kick mechanism (see Cordes 
& Chernoff 1998). For example, the magnetic field strengths required for 
the neutrino-driven mechanisms are >10 15 G, much larger than the currently 
inferred dipolar surface fields of typical radio pulsars; the internal magnetic 
fields of neutron stars and their evolution remain clouded in mystery; and 
several different mechanisms may contribute to the observed kick velocities. 

It is of interest to note that soft gamma repeaters ("magnetars": neutron 
stars with observed magnetic fields >10 14 G; see Thompson & Duncan 1996) 
have very high velocities, >1000 — 2000 km s _1 (e.g., Kaspi 1999; Marsden 
et al. 1999). Such a high velocity may well require supers trong magnetic fields 
(>10 16 G) to be present in the proto-neutron stars, although hydrodynamical 
effects remain a viable kick mechanism if enough presupernova asymmetry 
can be generated. It has recently been suggested that MHD jets can play an 
important role in supernovae (Khokhlov et al. 1999), but the origin of the jets is 
unknown, nor it is clear why the two opposite jets are so different (a necessary 
condition to produce a kick). 

I thank my collaborators Phil Arras, Peter Goldreich and Yong-Zhong Qian 
for their important contribution and insight. This work is supported by NASA 
grant NAG 5-8356 and by a fellowship from the Alfred P. Sloan foundation. 
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